Spin orientation and Spin-Hall effect induced by tunneling electrons 
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It is shown that a flux of unpolarized electrons across a symmetric double barrier quantum well 
induces a spin polarization inside the well. Besides, the transmitted current acquires a spin polarized 
component and the spin-Hall current flows in the planar direction. These phenomena are due to 
a combined effect of Dresselhaus interaction and the spin-orbit interaction induced by gradients of 
heterostructure material parameters. In contrast to previous studies of the spin filtering effect, we 
predict that it can be observed in case of an isotropic distribution of incident electrons. 
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I. INTRODUCTION 

The spin-orbit interaction (SOI) is a fundamental 
quantum relativistic phenomenon which recently at- 
tracted much interest in connection with spin transport 
of electrons in semiconductors and metals. Due to SOI, 
an electric field can influence the spin degree of freedom, 
thereby giving rise to a number of transport phenomena 
which have potential for application in spintronics. One 
of them is the spin-Hall effect (SHE) , which recently has 
been intensively studied both theoretically (for a review 
see[l|) and experimentally i^i^, A standard system to study 
this effect is a 2D electron gas confined within a quan- 
tum well. Due to SHE the electric current in the quan- 
tum well gives rise to a perpendicular flux of the spin 
polarization, as well as to the out of plane spin density 
near sample edges. Closely related to SHE is the electric 
spin orientation, which is a bulk inplane spin polariza- 
tion induced by the DC electric current.^ In both cases 
SOI is realized either via impurity scattering, or due to 
an intrinsic spin-orbit coupling mechanism. The latter 
consists of two parts. The first one is the Dresselhaus 
interaction^, which is inherent to all zinc-blende semicon- 
ductors. The second contribution is determined by gra- 
dients of material parameters and the electric potential 
across a heterointerfacei^ In quantum wells (QW) these 
interactions are averaged with wave functions of confine- 
ment coordinate in the heterostructure growth direction. 
After such an averaging procedure the second term trans- 
forms into the Rashba SOI® which is not zero only in 
asymmetric in z-direction (growth direction) heterostruc- 
tures. The Rashba and averaged Dresselhaus interactions 
are basic SOI widely used in works on SHE, as well as in 
works on spin-dependent transport in general. 

In this work we will consider spin-orbit effects on elec- 
tron transport from a different point of view. Namely, we 
will consider the electric current parallel to the growth 
direction, rather than parallel to the quantum well. An 
appropriate model for studying such a situation is a dou- 
ble barrier quantum well. A principal difference from 
the conventional 2D system is that one can not use nei- 
ther Rashba nor Dresselhaus interactions averaged over 



QW confinement. For example, a part of the latter in- 
teraction, which is proportional to the z-componcnt of 
the electron momentum operator, turns to when av- 
eraged with a wave function of a confined state in the 
well. At the same time it is finite for tunneling states. 
Considering these states it is also easy to see that SOI 
associated with electric fields at heterointerfaces does not 
reduce to the single parameter Rashba interaction. The 
explicit dependence on z of the initial spin-orbit interac- 
tion becomes important. To make this point more clear 
we deliberately considered a symmetric heterostructure, 
where Rashba SOI is zero. For such a model of a sym- 
metric double barrier quantum well we found out that 
the electric current in the z-direction induces a parallel 
to the z-axis spin density, as well as a current of spins 
polarized in a planar direction. The spin current flows 
parallel to heterointerfaces. The former effect is an ana- 
log of the electric spin orientation, while the latter is the 
spin-Hall effect. Besides, we found the transmittance of 
the double barrier structure to be dependent on the spin 
orientation. A similar spin filtering effect has been con- 
sidered before within various modelsi^i^ It was shown 
there that this effect can be observed only in case of 
an anisotropic distribution of incident electrons. It was 
suggested to create such a distribution applying a pla- 
nar electric field. In contrast to these works, we predict 
the spin filtering effect for an isotropic distribution of 
electrons tunneling through a symmetric double barrier 
structure. Moreover, our spin filter makes electrons to 
be polarized in z-direction, instead of the planar polar- 
ization in Ref. [Miol. Such fundamental distinctions arise 
from the ~ kz term of the Dresselhaus interaction which 
has been neglected in Ref. ITol . 

The general Hamiltonian of the problem will be derived 
in Sec. II. In Sec. HI we will present our results related to 
the spin-Hall effect, spin orientation and the spin filtering 
effect. A brief conclusion is presented in Sec. IV. 



II. SPIN-ORBIT HAMILTONIAN 

Let us consider a quantum well (QW) of the width 2d 
separated from the left (z < —d — h) and right (z > d + h) 
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parts of a doped semiconductor system by two equal bar- 
riers of the thickness b. These parts are assumed to be 
thermal reservoirs with respective chemical potentials ^/ 
and fir- Besides a system homogeneous in x,y directions, 
we will also consider an electron gas confined in the y- 
direction. For simplicity, the adiabatic case will be con- 
sidered when the confinement width w slowly increases 
from QW towards reservoirs. This situation is realized 
when the confinement is achieved by depleting the elec- 
tron gas with the help of electrodes on top of QW. 

The spin-orbit interaction in such a system is repre- 
sented by two Hamiltonians and Hso^ . The former 
is the Dresselhaus SOI, while the latter is SOI due to 
change of the band gap width and other material pa- 
rameters across heterointerfaces. Usually, in narrow gap 
semiconductors Hso is much stronger. Therefore, we will 
consider the Dresselhaus interaction within the first or- 
der perturbation theory. Only the part of -ffio^ which is 
proportional to the z-component of the electron momen- 
tum operator will be taken into account. Hence, for [001] 
growth direction 
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where (Tz is the Pauli matrix and 7 is the coupling param- 
eter, which is assumed to be z-indepcndent. An impor- 
tant property of Hamiltonian ([1]) is that its expectation 
values taken with tunnelling states incident from the left 
and from the right reservoirs have opposite signs, while 
it does not change sign when kx,ky —^ —kx,—ky. Alter- 
natively, other parts of the Dresselhaus SOI, which have 
been omitted in ([1]), are even functions of k^ and odd 
functions of kx,ky. Since, the interaction represented by 

h'^V is of the same symmetry, the omitted terms 

o 

do not add particularly new qualitative features to spin 
dependent electron tunneling, as well as to other effects 
considered below. At the same time, the symmetry dif- 
ference of ([1]) and Hso has important consequences for 
these effects. This is also the main reason why the results 
of Ref. [13, where interaction ([T]) has been neglected, are 
qualitatively different from those presented below. 

Following Ref. the Hamiltonian can be written 

as 



where 



■f^io^ = 7- {'^xky - (Jykx) h{z) , 



(2) 



(3) 



and fc|| = w kl 



Ky. 



The parameter h(z) denotes SOI 

strength which varies across the heterostructure depend- 
ing on semiconductor material parameters and the elec- 
tric potential. Ignoring the electron energy, which is 
much less than the gap value, /3(z) can be written as^ 

1 A(z) 



where Eg{z) and A(z) are respective values of the band 
gaps and split off energies. 

The major effect of SOI ^ is that it gives rise to spin 
precession around the z-axis. This precession takes place 
during particle transmission through the double barrier 
structure. Since the width of this structure is small in 
comparison with the spin precession length, the effect 
of the spin precession is expected to be small. In order 
to get explicitly the corresponding small parameter, the 
Hamiltonian can be transformed using an appropriate 
unitary transformation. Taking into account that the 
kinetic energy operator in z-direction is 
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2 dz m{z) dz 

one can apply the unitary transformation 



with 



where 



U 



(5) 



(6) 



(7) 



/3(z) 



2TO(z)3£;g(z)-|-A(z) 



(4) 
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This transformation removes ([1]) from the Hamiltonian. 
At the same time, applying it to ^ one obtains, up to 
the linear in 'd[z) terms, the spin orbit interaction 

Hso^H^f + ^^{axkx + ayky)h{z). (9) 

For the model under consideration, with rectangular 
symmetric barriers and a rectangular QW, h{z) becomes 

h{z) = k\\[{Pr- M{Siz-b^d)-6iz + b + d)) 

+ iPb- f3^)iSiz-d)-S{z + d))]. (10) 

The parameters Pr, Pb and /3„ denote SOI strengths for 
reservoirs, barriers and QW, respectively. 

The transmission wave functions are represented by 
two sets of functions incident from the left (-00 and from 
the right {ip^) of the double barrier structure. In the 
zeroth order, when the second term in ([9]) is ignored, 
these functions can be conveniently written using the chi- 
ral basis. In case of the quantum wire confinement this 
basis corresponds to the spin quantization axis directed 
along the y-axis. Outside the double barrier structure the 
scattering eigenstates are represented by incident, trans- 
mitted and reflected plane waves, with transmission and 
reflection amplitudes ta^ and rj^ , respectively, where 
cr = 1,2 denotes the spin projection in the correspond- 
ing chiral basis. This projection is conserved upon the 
scattering, as far as the second term in ^ is neglected. 
The wave vector of the scattering states is denoted as 
k — ^2'mj.{E — -Ell), where is the electron effective 
mass in the left and right reservoirs, E is the total energy 
and E\\ is the energy of motion in x, y directions. 
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III. SPIN CURRENT AND SPIN 
POLARIZATION 

A. Spin current 

The non zero spin current J*, where s and n denote 
the spin polarization and current direction, respectively, 
can be calculated already in the zeroth order with re- 
spect to "diz), while the latter will be important below in 
calculation of the spin density and spin dependent trans- 
mission. Hence, in this subsection we entirely neglect the 
presumably weak Dresselhaus SOI. Let us take the chi- 
ral component of J^^. That means that we look for the 
flux of spins polarized perpendicular to the flux direction, 
that can be expressed as = Je''"^, where e*"^ is the 
antisymmetric tensor. Using the conventional definition 
of the spin current operator — {w„,crs/4}, where the 
spin dependent part of 



m{z) 



h{z) 



(11) 



is obtained from ([2]) and ([3]) , the spin current density can 
be written as 



~snz /"oo 



dk , fcii 



k||,7 

fcii 



(12) 



where ^ — l,r and v?p[E) is the Fermi distribution func- 
tion for the left and right reservoirs. The first term in 
square brackets represents the bulk spin current density 
distributed in QW, barriers and outside, while the sec- 
ond term is the "surface" term which, according to pO)) . 
is finite only on heterostructure interfaces. From Eq. fT^ 
it becomes immediately evident that the spin current is 
not zero in the equilibrium state when np{E) — n'^p{E). 
For example, the surface current at each interface is given 
by p{z)AP(z), where p{z) is the local equilibrium elec- 
tron density and A/3(z) is the difference of the spin-orbit 
coupling parameters /3 on both sides of the interface, as 
follows from (fTO|) . In a symmetric QW the surface cur- 
rents on opposite interfaces flow in opposite directions, 
so that they cancel each other. It is easy to see that 
the total equilibrium current, obtained by integration of 
Eq. (fT2l) over z, is identically zero in case of a symmetric 
heterostructure. That follows from the symmetry rela- 
tion 'ij//^{z) = ip2^^{—z) which makes J^{z) to be an odd 
function of z. At the same time, one can not expect the 
total current to be zero in an asymmetric heterostruc- 
ture, as has been shown by Rashba^^ for confined states. 
It should be noted that Rashba found that the total equi- 
librium current is cubic with respect to the spin-orbit 
coupling constant, while the current density given by Eq. 
(IT^ is linear. The latter becomes evident from the above 
expression for the surface current. That means that, at 



least, linear and quadratic terms vanish after integration 
of dni) over z. 

By convention, the " nonequilibrium" current could be 
defined as a part of (|12p which is proportional to nlp — nTp. 
In symmetric QW this current density is an even function 
of z and, hence, the corresponding total nonequilibrium 
current is finite. However, one can not define unambigu- 
ously the dissipative part of the spin current using only 
its definition (|12p . Calculation of the spin accumulation 
at the sample boundary would be helpful to clarify the 
physical meaning of Eq. ([T! 



B. Spin orientation 

In order to calculate the spin density induced by the 
tunneling current, the second term in © must be taken 
into account. It causes spin flip processes upon transmis- 
sion and reflection of particles incident onto the double 
barrier structure. Therefore, we label spin variables of 
wave functions by two indices, as V'q/3j where a denotes 
the spin polarization of the incident wave. Treating such 
functions as matrices, the spin density can be expressed 
as 



1 dk 



(13) 



where cr = (ctj., (Ty,(7z) is the vector of Pauli matrices. We 
calculate (1131) in the first order perturbation theory with 
respect to i?(z). Since a commutator of two terms in ([9|) 
is proportional to <Tz , one should expect the z-component 
of S(z) to be finite. Further, the time inversion symmetry 
dictates 

E J^K^C'^^kj = - E Trir.iar.u,] ■ (w) 

Applying this relation to Eq. (fT^ the latter is trans- 
formed to 

poo Ih, 

^^(^) = E / j-Tr[i^'^az^P']{nUE)-nUE)). (15) 



It immediately follows from this expression that the spin 
density is zero in the equilibrium state. 

The first order correction to the wave function can be 
written in terms of the retarded Green function, so that 



i^Lpi^) = i^Liz)So.p + J dz'G0{z,z')Vpo.{z'Wc,{z'), 

(16) 

where Vfja{z') is a matrix element of the second term in 
^ and ipa{z) is the unperturbed wave function. In its 
turn, the Green function is given by 



Gpiz,z') = _^^[4(z)V>5(z')0(z-z') 
+ i;'p{z)^'^{z')9{z' - z)]. 



(17) 
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In the considered case of a symmetric heterostructure, 
it is easy to see that the transmission coefficient ti — 
t2 = t. At the same time, the reflection is spin depen- 
dent through its phase. We note that, according to ([H]), 
the matrix elements Vpa are proportional to k"^ — ky. Af- 
ter integration in (fT3|) over angles of the vector k|| this 
expression turns to 0. It does not necessary happen for 
other than [001] crystal orientations. We will not con- 
sider such an opportunity here. Also, we will not discuss 
here other than Dresselhaus's SOI effects, for example, 
due to strain, or due to potential gradients along x, y 
axes. Instead, let us consider a situation when electrons 
are confined, say, in y-direction. In this case k^ — ky be- 
comes kl — kl , where the overline and the label n denote 
averaging of the momentum operator over n-th quantum 
eigenstate in y-direction. Thus, the symmetry between x 
and y directions is broken and Sz does not turn to zero. 
Below we will assume a parabolic confinement. After av- 
eraging over the y-direction the spin-orbit interaction ([9|) 
depends only on kx, with fc|| — \kx\. 

Substituting ^ into ^ and then into Eq. ^ 
we express Sz{z) through unperturbed eigenstates. The 
latter are calculated for a square barrier structure de- 
scribed above. For an order of magnitude evaluation of 
the spin density, the result can be written in an ana- 
lytical form. Calculations are strongly simplified if only 
resonance terms in (|15|) are taken into account. One may 
also make use of the small parameter fc^,/2mu,J7, where 
kyj is the wave vector in the z-direction within QW and 
U is the barrier height. At the lowest transmission res- 
onance ku, — ko K n/2d. By this way, in the leading 
approximation the spin density in the center of the well 
(z=0) and the center of the wire (y=0) can be written as 



Sz =- -El'^"(o)p 



dk 



rrirmb 







[{ku 



-Hr2]2 Kk 



X A'&nid)h cosh Kbsmh2Kb 



(18) 



where h = kx{(3w — Pb-, k « y/2mbU. The phase 
•dn{z) is obtained from Eq. ([5]) by averaging ky with 
oscillatory wave functions (j)n{y)- The width of the 
transmission resonance in k space is given by F = 
{m^/mrmw){kkw / K^d) shih~'^ Kb^ where mr,mb,mi^ are 
effective masses in reservoirs, barriers and QW, respec- 
tively. A is SL dimensionless function of k. The value 
of this function is close to 1 in the range of parameters 
under consideration. 

For a numerical evaluation the following parame- 
ters have been taken: d=100A, 6=40A, electron den- 
sity in reservoirs 71=10^^ m~'^, the quantization energy 
of the parabolic confinement huj = A meV, 7=27 eV 
A'^F^ Other parameters correspond to the InAs quan- 
tum well, Ino,53Gao,47As barrier, and Ino.gGao^iAs reser- 
voirs. With such parameters we calculated F= 10''^ cm~^ 
and -did) w 10^''. From (|18p . in the linear regime 
A/i = — /i,. <ti Ep, the spin density can be evaluated 



as Sz ~ 0, lA/i meV~^ fim~^. Although the spin density 
is not high, nevertheless, in the range of 1 mV terminal 
voltages it can be detected by the Kerr rotation method^. 
The parameters of the heterostructure can also be opti- 
mized to reach the higher density. Experimentally, the 
total spin polarization can be enhanced in superlattices. 



C. Spin dependent transmittance 

The second term in ([9]) causes spin flips of a particle 
transmitting through the double barrier structure. The 
spin dependent transmission is obtained from (|17p and 
(fTH]) where tp^ is calculated a,t z > d + b. By this way, 
near the resonance, the spin flip transmittance becomes 



Atl,2 = 
At2A = 



-il6t^^{d)B , 
mi 



h'^kl, cosh^ nb 



K^k sinhnb 



(19) 



where S ~ 1 is a dimensionless factor. It is easy to see 
that, according to (|19p . Ati_2 transforms an unpolarized 
flux of electrons, say, from the left reservoir into a flux 
of spin polarized electrons, with the polarization in z- 
direction. Indeed, according to our choice of the spin 
basis, the z-polarization is obtained by taking an average 
of the ay spin operator. Eq. ifTO)) gives rise to two spinors 
corresponding to two possible initial spin polarizations: 
ipi — (i, Aii^2) and tjj2 — (Ai2,i,i)- For an unpolar- 
ized source, after summing up averages of ay with these 
spinors, one obtains a finite result, taking into account 
that near the resonance t = iT[kw — kg + iT]~^ w 1. 
It should be noted, that besides Aii_2 contains one 
more term which, however, does not result in a resonant 
polarized transmission. This term has been neglected. 
The transmission coefficient obtained from (jl9p is sym- 
metric with respect to x — x, y ~y and, hence, 
remains finite after the angular averaging, in contrast to 
a spin filtering effect considered in Ref 9,10. Close to the 
Fermi energy and with the same numerical parameters 
used for evaluation of (|18p we obtain Aii^2 ~ 3 • 10~^. 
In spite of its small value, this coefficient can lead to 
a noticeable accumulation of the spin polarization in a 
reservoir with small spin relaxation rate. For example, 
let us take a reservoir of 1 fim!^ volume. The spin po- 
larization flux carried from the left reservoir through the 
wire of the length L is given by 



E 



dk—lm[t*Ati,2][np{E) 



Taking iii — fir = ImeV, L = 1/im, the typical spin re- 
laxation time in bulk semiconductors 1 ns and all the 
rest of the parameters as in above evaluations, we get the 
spin density in the reservoir around 0,5 fim~^, which is 
within the optical detection rangei^ 
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IV. CONCLUSION 

We considered the spin orbit effects associated with 
resonant tunnehng of electrons through a symmetric dou- 
ble barrier structure. Two contributions to SOI have 
been taken into account, namely, the Dresselhaus SOI 
and the spin-orbit interaction induced by gradients of 
heterostructure material parameters. We found out that 
the vertical transport of electrons gives rise to the spin 
current flowing parallel to heterointerfaces, as well as to 
the spin polarization within QW. These effects are anal- 
ogous to the spin-Hall effect and the electric spin orienta- 
tion, intensively studied recently for 2DEG. A distinction 
with these traditional studies is that instead of electron 
wave functions confined in QW we employ eigenstates 
corresponding to resonant electron transmission. More- 
over, instead of an electric field parallel to QW, the ver- 
tical bias has been considered. In compliance with such 
a 3D model, we calculated a distribution in z-direction of 
the spin-Hall current density and spin polarization. This 
dependence reveals interesting structure, such as "sur- 
face" currents flowing along the heterointerfaces. The 
spin-Hall current density does not turn to with the 
external bias, thus signalling about existence of the equi- 
librium spin current density, which was found to be linear 
with respect to the spin-orbit interaction. At the same 



time, in a symmetric double well structure the net equi- 
librium current, obtained by integration of the current 
density over z, turns to 0. Although the spin-Hall and the 
equilibrium currents take place in absence of the Dressel- 
haus SOI, the latter is necessary to obtain the finite out of 
plane spin density within QW. It was found out that the 
most important is a part of the Dresselhaus SOI which is 
proportional to the z-component of the electron momen- 
tum operator. This SOI also gives rise to spin dependent 
transmission. Due to interplay of the Dresselhaus SOI 
and the spin-orbit interaction induced by gradients of 
heterostructure material parameters, an unpolarized and 
isotropic in x, y directions beam of electrons becomes po- 
larized in the z-direction after tunneling through double 
barrier QW. 

The effects discussed in this paper can be interesting 
in application to metal surfaces with strong spin-orbit ef- 
fects associated with surface states. The STM set up 
with a magnetic tip could be employed for studying spin 
related effects. The above model, however, must be mod- 
ified to take into account spin-orbit interactions typical 
for a particular metal surface. 
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